Gapless Fermi Surfaces in anisotropic multiband superconductors in magnetic field. 
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We propose that a new state with a fully gapless Fermi surface appears in quasi-2D multiband 
superconductors in magnetic field applied parallel to the plane. It is characterized by a paramagnetic 
moment caused by a finite density of states on the open Fermi surface. We calculate thermodynamic 
and magnetic properties of the gapless state for both s-wave and d-wave cases, and discuss the details 
of the 1-st order metamagnetic phase transition that accompanies the appearance of the new phase 
in s-wave superconductors. We suggest possible experiments to detect this state both in the s-wave 
(2-H NbSe2) and d-wave (CeCoIns) superconductors. 
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In the original BCS theory of superconductivity (SC) 
excited states are separated from the ground state by 
an energy gap. SC does not necessarily lead to a fully 
gapped energy spectrum of quasiparticle excitations. It 
is well known that for unconventional SC, where the 
mechanism of SC is different from the BCS, non-trivial 
symmetry of the order parameter allows the existence of 
point or line nodes on the Fermi surface (FS) in momen- 
tum spacep], |H, H|. For ordinary s-wave pairing gapless 
SC appears in the presence of paramagnetic impurities, 
when the time- reversal t — > — t symmetry is broken^. 
Neglecting the orbital effects, for a SC placed in mag- 
netic field one expects the realization of inhomogeneous 
Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) state [a, [1], in 
which the superconducting gap also passes through zeros 
this time at the points in real space. 

A resurgence of interest to the new possibilities of gap- 
less SC is related to the recent proposal of "interior gap 
superfluidity" of Liu and Wilczek 7] , in which whole re- 
gions of the Fermi surface remain ungapped. Liu and 
WilczekQ have considered a situation often realized in 
atomic physics and high-energy physics fl[, and the Bose- 
Einstein condensation, in which the superfluid pairing 
takes place for non-identical fermions condensed by an 
optical trap. In this case (often called unbalanced (UB) 
pairing), due to both the difference in concentrations and 
bare masses, the two condensing components would have 
to form two FS with different Fermi momenta. For SC 
pairing between the two different Fermi surfaces it is en- 
ergetically favorable to leave a part of momentum space 
near the larger FS in normal state. Ground states of 
this type have been studied for ordinary SC some time 
ago[iJ[l(|. However, the tendency toward inhomogeneous 
LOFF state in ordinary s-wave SC always prevailed. 

In this letter we show that similar features in the 



energy spectrum should naturally appear in quasi-2D 
multiband SC, such as some organic SC or "115" heavy 
fermion materials, CeMIns (M = Co,Ir). Most recent 
experimental activity has been devoted to finding the 
LOFF state [H E 111 0, El in CeCoIn 5 . The latter 
material is characterized by quasi- 2D Fermi surfaces and 
a multiband energy spectrum. We consider a model of 
quasi-2D two-band SC of both s-wave and d-wave type in 
magnetic field applied parallel to the plane. Apart from 
the LOFF state in higher magnetic fields, for hbB = I 
comparable with the smaller gap we observe whole re- 
gions of open FS, similar to the situation considered in 
Ref. 0. For an s-wave SC we investigate analytically in 
detail this low-temperature and low field region of the 
phase diagram. 

We adopt the standard multiband interaction scheme 
(see, e.g., Ref.[16]). The matrix elements Uik{p;p') for 
the interaction enter the definitions of the gaps, Aj(p) 
for each Fermi surface (FS) as: 

Ai(p) = -Tj2 U lk (p,p')F k (iLu n ,p'), (1) 

n,k,p' 

where Fk(iu> n ,p) is the anomalous Gor'kov function, k,i 
are band indices, Uik is the interaction between bands i 
and fc. 

The type of the superconducting state below T c de- 
pends on the choice of the pairing ansatz : 



Uik(p;p') = xivWikxiv') 



(2) 



where x(v) i s the appropriate irreducible representation; 
we take x(</?) below as a const (1) for the s-wave pair- 
ing, or as cos (2ip) for the d-wave pairing. Solutions of 
the multiband Gor'kov equations \yA Il8j for the Green's 
functions in magnetic field I = fJ-sB can be written as: 
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(3) 
(4) 
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The energy spectrum of the system for excitations near 
each FS is given by the poles of the Gfc(u;„,p): 



E k (p) = V6(P) 2 + |A fe (p)| 2 + la* 



(5) 



The bands with different k are coupled by the gap equa- 
tion, Eq.(TT]). We consider below a model with 2 FS. 

While the search for the LOFF-state commonly starts 
from the side of higher fields, we study the effects in small 
magnetic fields of the order of the smaller gap, A2. Our 
main interest lies in the field range where A2 < / < Ai. 
Once the magnetic field I = fisB exceeds the smaller 
gap, then, according to Eq. (f5|) , electron- and hole- pock- 
ets will open, forming an ungapped area near the second 
Fermi surface (FS2). In the s-wave case this process is 
accompanied by a weak 1st order phase transition. For 
a d-wave SC the energy gaps Afc(p) have line nodes and 
associated gapless states from the start. Nevertheless, 
depending on the strength of interactions, an irregular 
behavior of the gap amplitude as a function of magnetic 
field also occurs in some region of model parameters (see 
Fig© > which indicates a 1-st order transition. 

It is convenient fl9j to express the solution of Eq.fTJ in 
terms of dimensionless coupling constants, A;fc = UikVk, 
where v% is the density of states on the k-th Fermi surface 
(FSk). The linearized gap equation Eq.JT]) leads to the 
familiarly instability curve for T c , which, we find, is 
independent of the number of FS involved: 



In — = * ( - I /?( 

T c0 \2 



T I 1 ■ 1 

2 2-KT r 



(6) 



Here T c q is the superconducting transition temperature 
without the magnetic field, 



T, 



cO 



(ff<0) 



A11A22 — A12A21 



An + A22 + v (An — A22) 2 + 4A12A 



(7) 
, (8) 
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7 ~ 1.781, and A is the upper cut-off for the interac- 
tions in Eq.([2|). However, the total (T,B) phase dia- 
gram for two bands changes significantly, especially at 
lower temperatures and fields. Some main qualitative 
changes in the physics of multiband SC in this area can 
already be seen in a simplified model with Ai as the pri- 
mar y g ap, and A2 <C Ai induced by the SC order on 
FSl[2l|. When I is close to the primary energy gap, Ai, 
an inhomogeneous LOFF state will appear [221 [231 ] . For 
Ai ~ A2 there could be significant modifications for the 
LOFF state, and the boundaries of the LOFF state on 
the (T, B) phase diagram. We assume that two gaps dif- 
fer enough for the LOFF state not to change significantly 
from the single-band model [13, [HI]. Below we consider 
in more detail the low-field region of the (T, _B)-plane. 

In the weak coupling approach, j « 1, the ratio of 
A 2 (T,B), the driven gap, and Ai(T,B), the primary 
gap, which we define as model parameter t, sets in at 
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FIG. 1: Magnetic field dependence of the primary energy gap 
Ai for the s-wave case for t = A2/A1 = 0.2, and different 
values of parameter a, 0.1, 0.2, 0.3, and 0.4. The energy gap 
at / = increases with growing a, see Eas(|9lTT2]l . The inset 
shows the corresponding metamagnetic transition. The calcu- 
lated paramagnetic moment is normalized to the full value of 
the normal moment from the second Fermi surface. The grey 
lines represent the positions of the 1-st order phase transition. 



(9) 



T c and is temperature- and magnetic field-independent 

A2 = 2Xi2 

Ai A22 — An + \J (An — A22) 2 + 4A12A21 

Eq. {JJ for the s-wave case can be easily solved analyt- 
ically at T = 0. Introducing new parameters, 

a = t 2 v 2 (v x + u 2 t 2 )-\ A = {n/j)T c0 , (10) 

we find two different solutions for Eq.((T]) for I < I cr — 
A (l + \/l - t 2 y a , and no solutions for / > I cr . The 
first solution is 



Ai = Am =t- a A Q , I< A 20 

— ( Ai V 

A \I+y/P-t 2 Al ) ' 



(11) 

A 20 < I < lev (12) 



Here A l0 = A,(T = 0,B = 0), % = 1,2. The second 
solution exists for magnetic fields Ao/2 < I < I cr , 

(I + ^I 2 - A\f- a {I + yfp - t 2 A\) a = A , (13) 

and is the familiar 0, [11] unstable solution for the energy 
gap in high magnetic fields. The two solutions are plotted 
in FigfTJ The re-entrant behavior in magnetic fields / ~ 
A20 = iAio clearly indicates the 1-st order character of 
transition into the gapless state. At this transition an 
open Fermi surface is formed, according to the energy 
spectrum given by Eq. ([5]) . The position of the 1-st order 
phase transition is found from the energy at T = 0, which 
for / < A20 has the usual form, 

AE = E s -E N0 = -{v 1 Aj + v 2 A 2 ,)/4. (14) 
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FIG. 2: Magnetic field dependence of the primary energy gap 
amplitude, Ai, for the d-wave case for t — A2/A1 = 0.2, and 
different values of parameter a, 0.1, 0.2, 0.3, and 0.4. The 
energy gap at I — increases with growing a. The plot is in 
reduced units of Aq = 2irT c o/"/^/e, the energy gap value cor- 
responding to T c o for a 1-band d-wave superconductor. The 
inset shows the corresponding metamagnetic transition, which 
only occurs at a > 0.24. The calculated paramagnetic mo- 
ment is normalized to the full value of the normal moment 
from the second Fermi surface. The grey lines represent the 
positions of the 1-st order phase transition. 

For / > A20, we also find a contribution from the normal 
excitations: 

AE = -(1/1A? + 1/2 A| + 2v 2 IiJp - Af)/4. (15) 

In Fig Q] we also show the field dependence of the para- 
magnetic contribution to the total magnetization for the 
same values of the parameter a. There is a characteris- 
tic metamagnetic jump in the magnetization at / = A20- 
At / > A20 the finite density of states (DOS) arises on 
each new electron- and hole- FS that will result in the 
linear- in- T specific heat at low temperatures. The 1-st 
order transition separating the two field regimes should 
be clearly seen in thermodynamic measurements. 

We have obtained similar results for d-wave multiband 
superconductors, such as the "115" materials. The gen- 
eral theoretical formulas involve an average over the an- 
gular variable tp, and are more cumbersome than Eqs 
(fTTj) - (|15p for the s-wave case. In Fig[2]we show, for com- 
parison, the dependence of the amplitude of the order pa- 
rameter and the metamagnetic transition for the d-wave 
pairing. Due to the presence of line nodes in d-wave case, 
the energy spectrum is gapless already at / = 0. Never- 
theless, similar processes take place near FS2 that could 
lead to a 1st order metamagnetic transition for some val- 
ues of the coupling constants. However, while for the 
s-wave case the 1-st order transition is always present, 
for d-wave it turns into a smooth crossover for a«l. 

The metamagnetic transition can be studied analyti- 
cally for the s-wave case when the second gap and the 



parameter a are small (2lT|. A20 <§C A10, a < 1. Intro- 
ducing 77 = (I - A 2 o)/A 2 o and r A ee (A 2 - A 2 o)/A 2 o, 
we find that in the vicinity of this transition Eq. (fT2]) is 
considerably simplified, 

n = (l/2)r| a - 2 +r A . (16) 

Expanding the energy Eq. (|15[) in the vicinity of this tran- 
sition, we find: 

E s - E NQ = E (l - 2r 2 - (4/3)r A a- 2 ), (17) 

where Eq is the energy of superconducting state at / = 0. 
The cubic terms in the energy, and the form of Eq (|16[) 
clearly indicate a first order transition. After a simple 
calculation, we find that the 1-st order transition occurs 
at Ti cr = — 3a 2 /8. The energy gap t a changes abruptly 
from r A = to t a = — 3a 2 /2, which corresponds to 
a metamagnetic transition, with a paramagnetic jump 
in the magnetic moment M = (3a/ 2) (13^2^.20 (or to a 
sudden appearance of the finite DOS). We find a simple 
expression for the magnetic moment for / > I cr in the 
vicinity of this transition: 

M = ^ B v 2 A 20 {a + yja 2 + 2tj). (18) 

The paramagnetic moment in the s-wave case always 
appears at / > A20 in a 1-st order phase transition. Nev- 
ertheless, in case of the driven second gap, A20 <C A10, 
the 1-st order transition is weak, and so is the correspond- 
ing change in the SC order parameter, A2 = tA±. In 
the first approximation this change can be neglected [2l|. 
Then the temperature- and field- dependence of the mag- 
netic moment is completely described by the standard [2J| 
formulas that follow from the energy spectrum Eq.([5]), 
where A2 = A20 is regarded as a constant. For example, 
a simple analytic expression for the magnetic moment in 
the s-wave case at T = is: 

M = n B v 2 y/l 2 - A 2 ,, I>A 2 . (19) 

Note that M = for I < A 2 . For a d-wave case, FS2 
gives a contribution for I above and below A2: 

M 2 = -[i B v 2 [ dtpJp-Al sinV, (20) 
t Jo v 

where the upper limit is A = n/2 for / > A2, or A — 
arcsin(I/A 2 ) for I < A 2 . This is an elliptic integral 
of second kind. Note that in the d-wave case FS1 gives 
the usual nodal contribution, Mi — 0.5/is^i/ 2 /Ai. The 
density of states for FS2 in the s-wave case is also given 
by a simple formula: 

y 2 (I) - u 2 I/^P -Al (21) 

In the d-wave case one has to introduce the appropriate 
angular averages of this result, and add the familiar nodal 
contribution from the first Fermi surface [24|. 
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In summary, we have shown that a gapless Fermi 
spectrum characterized by open Fermi surfaces is an in- 
evitable feature for a quasi-2D multiband superconductor 
placed into a large enough field parallel to the plane. The 
new state is fully analogous to the one studied in Ref.Q 
for the unbalanced pairing problem. Unlike Ref. [3] , how- 
ever, such a gapless state sets in as the 1-st order tran- 
sition in increased magnetic field. Measurements of the 
specific heat in applied field are the most direct way to 
observe the effect in s-wave superconductors, such as 2H- 
NbSe2 [HI, [2(| . The transition also leads to a metamag- 
netic jump in the magnetization. For a d- wave pairing, 
because of the nodes, gapless excitations are present even 
without external field. As the field is increased, the open 
Fermi surfaces develop gradually, although character of 
the process may depend on the interaction parameters. 
Applied fields should be low enough for these phenom- 
ena not to interact with the LOFF state. It is broadly 
believed [111, [T3, EH that the properties of CeCoIns may 
be close enough to a two-dimensional model to display 



the inhomogeneous LOFF state [27|- If CeCoIns, indeed, 
belongs to the strongly quasi-2D class, the low field prop- 
erties studied above should manifest themselves as well. 
Then, if interactions in CeCoIns were strong enough to 
result in a 1-st order transition of Fig[2l the latter could 
be observed best by calorimetric measurements, as for the 
s-wave pairing. If not, then one may rely on the NMR 
methods for the observation of a rather non-monotonic 
field behavior for non-linear susceptibility shown in the 
insert of Fi^5] (we have not considered possible implica- 
tions of the effect for thermal conductivity in the presence 
of a magnetic field) . The above effects should be expected 
for other 2D organic compounds. An ideal realization of 
the scheme would be superconductivity localized at the 
surface . 
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